Atmospheres and wind properties of non-spherical AGB stars

The wind-driving mechanism of asymptotic giant branch (AGB) stars is commonly attributed to a two-step process: first, gas in the stellar atmosphere is levitated by shockwaves caused by stellar pulsation, then accelerated outwards by radiative pressure on newly formed dust, inducing a wind. Dynamical modelling of such winds usually assumes a spherically symmetric star. We explore the potential consequences of complex stellar surface structures, as predicted by three-dimensional (3D) star-in-a-box modelling of M-type AGB stars, on the resulting wind properties with the aim to improve the current wind models. Two different modelling approaches are used; the CO$^5$BOLD 3D star-in-a-box code to simulate the convective, pulsating interior and lower atmosphere of the star, and the DARWIN one-dimensional (1D) code to describe the dynamical atmosphere where the wind is accelerated. The gas dynamics of the inner atmosphere region at distances of $R\sim1-2R_\star$, which both modelling approaches simulate, are compared. Dynamical properties and luminosity variations derived from CO$^5$BOLD interior models are used as input for the inner boundary in DARWIN wind models in order to emulate the effects of giant convection cells and pulsation, and explore their influence on the dynamical properties. The CO$^5$BOLD models are inherently anisotropic, with non-uniform shock fronts and varying luminosity amplitudes, in contrast to the spherically symmetrical DARWIN wind models. DARWIN wind models with CO$^5$BOLD-derived inner boundary conditions produced wind velocities and mass-loss rates comparable to the standard DARWIN models, however the winds show large density variations on time-scales of 10-20 years.Comment: 13 pages, 12 figures, Accepted for publication in A&

Pulsation-induced atmospheric dynamics in M-type AGB stars. Effects on wind properties, photometric variations and near-IR CO line profiles

Wind-driving in asymptotic giant branch (AGB) stars is commonly attributed to a two-step process. First, matter in the stellar atmosphere is levitated by shock waves, induced by stellar pulsation, and second, this matter is accelerated by radiation pressure on dust, resulting in a wind. In dynamical atmosphere and wind models the effects of the stellar pulsation are often simulated by a simplistic prescription at the inner boundary. We test a sample of dynamical models for M-type AGB stars, for which we kept the stellar parameters fixed to values characteristic of a typical Mira variable but varied the inner boundary condition. The aim was to evaluate the effect on the resulting atmosphere structure and wind properties. The results of the models are compared to observed mass-loss rates and wind velocities, photometry, and radial velocity curves, and to results from 1D radial pulsation models. Dynamical atmosphere models are calculated, using the DARWIN code for different combinations of photospheric velocities and luminosity variations. The inner boundary is changed by introducing an offset between maximum expansion of the stellar surface and the luminosity and/or by using an asymmetric shape for the luminosity variation. Models that resulted in realistic wind velocities and mass-loss rates, when compared to observations, also produced realistic photometric variations. For the models to also reproduce the characteristic radial velocity curve present in Mira stars (derived from CO $\Delta v = 3$ lines), an overall phase shift of 0.2 between the maxima of the luminosity and radial variation had to be introduced. We find that a group of models with different boundary conditions (29 models, including the model with standard boundary conditions) results in realistic velocities and mass-loss rates, and in photometric variations

Exploring wind-driving dust species in cool luminous giants III. Wind models for M-type AGB stars: dynamic and photometric properties

Stellar winds observed in asymptotic giant branch (AGB) stars are usually attributed to a combination of stellar pulsations and radiation pressure on dust. Shock waves triggered by pulsations propagate through the atmosphere, compressing the gas and lifting it to cooler regions, which create favourable conditions for grain growth. If sufficient radiative acceleration is exerted on the newly formed grains through absorption or scattering of stellar photons, an outflow can be triggered. Strong candidates for wind-driving dust species in M-type AGB stars are magnesium silicates (Mg$_2$SiO$_4$ and MgSiO$_3$). Such grains can form close to the stellar surface, they consist of abundant materials and, if they grow to sizes comparable to the wavelength of the stellar flux maximum, they experience strong acceleration by photon scattering. We use a frequency-dependent radiation-hydrodynamics code with a detailed description for the growth of Mg$_2$SiO$_4$ grains to calculate the first extensive set of time-dependent wind models for M-type AGB stars. The resulting wind properties, visual and near-IR photometry and mid-IR spectra are compared with observations.We show that the models can produce outflows for a wide range of stellar parameters. We also demonstrate that they reproduce observed mass-loss rates and wind velocities, as well as visual and near-IR photometry. However, the current models do not show the characteristic silicate features at 10 and 18 $\mu$m as a result of the cool temperature of Mg$_2$SiO$_4$ grains in the wind. Including a small amount of Fe in the grains further out in the circumstellar envelope will increase the grain temperature and result in pronounced silicate features, without significantly affecting the photometry in the visual and near-IR wavelength regions.Comment: 11 pages, 14 figure

Tomography of silicate dust around M-type AGB stars I. Diagnostics based on dynamical models

The heavy mass loss observed in evolved asymptotic giant branch stars is usually attributed to a two-step process: atmospheric levitation by pulsation-induced shock waves, followed by radiative acceleration of newly formed dust grains. Detailed wind models suggest that the outflows of M-type AGB stars may be triggered by photon scattering on Fe-free silicates with grain sizes of about 0.1 - 1 $\mu$m. Due to the low grain temperature, these Fe-free silicates can condense close to the star, but they do not produce the characteristic mid-IR features that are often observed in M-type AGB stars. However, it is probable that the silicate grains are gradually enriched with Fe as they move away from the star, to a degree where the grain temperature stays below the sublimation temperature, but is high enough to produce emission features. We investigate whether differences in grain temperature in the inner wind region, which are related to changes in the grain composition, can be detected with current interferometric techniques, in order to put constraints on the wind mechanism. To investigate this we use radial structures of the atmosphere and wind of an M-type AGB star, produced with the 1D radiation-hydrodynamical code DARWIN. The spectral energy distribution is found to be a poor indicator of different temperature profiles and therefore is not a good tool for distinguishing different scenarios of changing grain composition. However, spatially resolved interferometric observations have promising potential. They show signatures even for Fe-free silicates (found at 2-3 stellar radii), in contrast to the spectral energy distribution. Observations with baselines that probe spatial scales of about 4 stellar radii and beyond are suitable for tracing changes in grain composition, since this is where effects of Fe enrichment should be found.Comment: Accepted for publication in Section 8. Stellar atmospheres of Astronomy and Astrophysics. The official date of acceptance is 07/09/2017. 9 pages, 7 figures, 4 figures in appendi

Dust-driven winds of AGB stars: The critical interplay of atmospheric shocks and luminosity variations

Winds of AGB stars are thought to be driven by a combination of pulsation-induced shock waves and radiation pressure on dust. In dynamic atmosphere and wind models, the stellar pulsation is often simulated by prescribing a simple sinusoidal variation in velocity and luminosity at the inner boundary of the model atmosphere. We experiment with different forms of the luminosity variation in order to assess the effects on the wind velocity and mass-loss rate, when progressing from the simple sinusoidal recipe towards more realistic descriptions. Using state-of-the-art dynamical models of C-rich AGB stars, a range of different asymmetric shapes of the luminosity variation and a range of phase shifts of the luminosity variation relative to the radial variation are tested. These tests are performed on two stellar atmosphere models. The first model has dust condensation and, as a consequence, a stellar wind is triggered, while the second model lacks both dust and wind. The first model with dust and stellar wind is very sensitive to moderate changes in the luminosity variation. There is a complex relationship between the luminosity minimum, and dust condensation: changing the phase corresponding to minimum luminosity can either increase or decrease mass-loss rate and wind velocity. The luminosity maximum dominates the radiative pressure on the dust, which in turn, is important for driving the wind. These effects of changed luminosity variation are coupled with the dust formation. In contrast there is very little change to the structure of the model without dust. Changing the luminosity variation, both by introducing a phase shift and by modifying the shape, influences wind velocity and the mass-loss rate. To improve wind models it would probably be desirable to extract boundary conditions from 3D dynamical interior models or stellar pulsation models.Comment: 11 pages, 13 figures, accepted for publication in A&

Three-component modeling of C-rich AGB star winds I. Method and first results

Radiative acceleration of newly-formed dust grains and transfer of momentum from the dust to the gas plays an important role for driving winds of AGB stars. Therefore a detailed description of the interaction of gas and dust is a prerequisite for realistic models of such winds. In this paper we present the method and first results of a three-component time-dependent model of dust-driven AGB star winds. With the model we plan to study the role and effects of the gas-dust interaction on the mass loss and wind formation. The wind model includes separate conservation laws for each of the three components of gas, dust and the radiation field and is developed from an existing model which assumes position coupling between the gas and the dust. As a new feature we introduce a separate equation of motion for the dust component in order to fully separate the dust phase from the gas phase. The transfer of mass, energy and momentum between the phases is treated by interaction terms. We also carry out a detailed study of the physical form and influence of the momentum transfer term (the drag force) and three approximations to it. In the present study we are interested mainly in the effect of the new treatment of the dust velocity on dust-induced instabilities in the wind. As we want to study the consequences of the additional freedom of the dust velocity on the model we calculate winds both with and without the separate dust equation of motion. The wind models are calculated for several sets of stellar parameters. We find that there is a higher threshold in the carbon/oxygen abundance ratio at which winds form in the new model. The winds of the new models, which include drift, differ from the previously stationary winds, and the winds with the lowest mass loss rates no longer form.Comment: 15 pages, 5 figures, accepted by A&

An extensive grid of DARWIN models for M-type AGB stars I. Mass-loss rates and other properties of dust-driven winds

The purpose of this work is to present an extensive grid of dynamical atmosphere and wind models for M-type AGB stars, covering a wide range of relevant stellar parameters. We used the DARWIN code, which includes frequency-dependent radiation-hydrodynamics and a time-dependent description of dust condensation and evaporation, to simulate the dynamical atmosphere. The wind-driving mechanism is photon scattering on submicron-sized Mg$_2$SiO$_4$ grains. The grid consists of $\sim4000$ models, with luminosities from $L_\star=890\,{\mathrm{L}}_\odot$ to $L_\star=40000\,{\mathrm{L}}_\odot$ and effective temperatures from 2200K to 3400K. For the first time different current stellar masses are explored with M-type DARWIN models, ranging from 0.75M$_\odot$ to 3M$_\odot$. The modelling results are radial atmospheric structures, dynamical properties such as mass-loss rates and wind velocities, and dust properties (e.g. grain sizes, dust-to-gas ratios, and degree of condensed Si). We find that the mass-loss rates of the models correlate strongly with luminosity. They also correlate with the ratio $L_*/M_*$: increasing $L_*/M_*$ by an order of magnitude increases the mass-loss rates by about three orders of magnitude, which may naturally create a superwind regime in evolution models. There is, however, no discernible trend of mass-loss rate with effective temperature, in contrast to what is found for C-type AGB stars. We also find that the mass-loss rates level off at luminosities higher than $\sim14000\,{\mathrm{L}}_\odot$, and consequently at pulsation periods longer than $\sim800$ days. The final grain radii range from 0.25 micron to 0.6 micron. The amount of condensed Si is typically between 10% and 40%, with gas-to-dust mass ratios between 500 and 4000.Comment: Accepted to A&A, 17 pages, 15 figure

Exploring wind-driving dust species in cool luminous giants II. Constraints from photometry of M-type AGB stars

The heavy mass loss observed in evolved asymptotic giant branch (AGB) stars is usually attributed to a two-stage process: atmospheric levitation by pulsation-induced shock waves, followed by radiative acceleration of newly formed dust grains. The dust transfers momentum to the surrounding gas through collisions and thereby triggers a general outflow. Radiation-hydrodynamical models of M-type AGB stars suggest that these winds can be driven by photon scattering -- in contrast to absorption -- on Fe-free silicate grains of sizes 0.1--1\,$\mu$m. In this paper we study photometric constraints for wind-driving dust species in M-type AGB stars, as part of an ongoing effort to identify likely candidates among the grain materials observed in circumstellar envelopes. To investigate the scenario of stellar winds driven by photon scattering on dust, and to explore how different optical and chemical properties of wind-driving dust species affect photometry we focus on two sets of dynamical models atmospheres: (i) models using a detailed description for the growth of Mg$_2$SiO$_4$ grains, taking into account both scattering and absorption cross-sections when calculating the radiative acceleration, and (ii) models using a parameterized dust description, constructed to represent different chemical and optical dust properties. By comparing synthetic photometry from these two sets of models to observations of M-type AGB stars we can provide constraints on the properties of wind-driving dust species. Photometry from wind models with a detailed description for the growth of Mg$_2$SiO$_4$ grains reproduces well both the values and the time-dependent behavior of observations of M-type AGB stars, providing further support for the scenario of winds driven by photon scattering on dust.Comment: Accepted for publication in A&A. 15 pages, 14 figure